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ABSTRACT 


ace bcemesOoun (i vicumm and Beaky. J.eAcoust. Soc. Am. v. 66, 1124-1130 
(1989)) that spilling breakers in the laboratory produce a sound close to the Knudsen 
sea surface noise spectrum from 400 ]]z to 20 klIz. The surface spectral production 
density of newly-created near surface bubbles under these laboratory spilling breakers 
has now been acoustically determined bv using an array of hvdrophones. The surface 
spectral density has been obtained by identifying the individual bubbles that create this 
noise spectrum. The radi were calculated from the resonance frequencies. The bubble 
positions on the surface of the water were determined from the difference in time of ar- 
rival of the bubble radiation to the elements of two vertical arrays of hydrophones. The 
production area and rate of production of bubbles of radi 0.048 to 7.40 millimeters have 
been calculated and the total volume of air encapsulated into bubble foam per unit area 


of spilling breakers has been determined. 
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I INTRODUCTION 


The purpose of this thesis is to develop a method for the determination of the bubble 
production density from breaking waves. The motivation for this study came from the 
observation that the spectrum produced from the noise of the breaking waves in the 
laboratory tank at the Naval Postgraduate School closely resembled the Knudsen sea 
surface spectrum of the ambient noise at sea. This development suggested that the 
breaking waves experimentally produced will exhibit the same bubble characteristics as 
breaking waves at sea. It 1s hoped that the capability of comparing the calculated bub- 
ble production density with the area of foam coverage monitored with satellite imaging 
and aerial photography will aid in acoustical and meteorological research of air/sea 
interactions by predicting the number and size of bubbles which mav be found in the 
region being studied. 

The method used in this research is to utilize passive acoustics to identify individual 
bubbles created from breaking waves. The breaking of the wave, characterized by the 
initual formation of a large amplitude. low frequency bubble, was used as a triggering 
source for the signal processing equipment. The location of each subsequent bubble was 
Saleulited by using the difference 1 tine of arrival of the noise from the bubble to the 
different elements of a hydrophone array. The bubble production area can be measured 
by positionmg all the bubble locations on a grid. The rate of bubble production was 
observed to behave as an exponentially decreasing function of time. This was determined 
by noting the relative time of arrival of the pressure waves from the individual oscillating 
bubbles to a hydrophone array. The bubble production density and volume rate of gas 
entrainment were calculated from the bubble population data and the measured area of 
bubble production. This technique provided information on bubble production on the 
surface, a quanuty which had not been directly measured in anv other previous exper- 


mental work. 


HW. BACKGROUND 


Many processes of acoustical, meterological, geochemical, and biological interest are 
associated with the presence of bubbles in the near surface of the oceans. Among these 
are ambient noise, aerosol formation, surface coagulation, fractionalization of organic 
or inorganic materials, and gas exchange. The bubble populations that influence these 
phenomena are composed of bubble patches produced by breaking waves superimposed 
on a background population of bubbles. 

The background population of bubbles originates from bubble patches due to pre- 
vious breaking waves. biological activity, droplet impacts, and sediment outgassing. 
Some of these bubbles will rise and eventually burst on the surface affecting the proc- 
esses already mentioned. Other bubbles undergo diffusion and will completely dissolve. 
An expression given by Levich (1962) for the rate at which these bubbles willidissmiie 


is described by the following equation: 


4 
at = Sasa aac 
where 1) 1s the number of moles of the gas. D 1s the gas diffusivity in wr] sec. U Ise 
bubble velocity in mi sec, a 1s the bubble radius m meters, and C 1s the Peclet mimes 
Wiaiclydesciibes ditusive aiass traiisier. 

The vertical distribution of the bubbles beneath the surface has been stidieaiins 
Thorpe (1986) and Farmer and Vagle (1989), Clouds of subsurface bubbles caused by 
wind waves breaking in deep water and creating whitecaps have been observed using a 
vertically pointing narrow beam sonar and a dual beam side seanning sonar. The vertical 
distribution of bubbles 1s influenced by several actions: the turbulent transport from a 
source at the surface. the rise of the bubbles due to buovancv. the change i bubBlemae 
dius due to the diffusion discussed earlier, and change of hydrostatic pressure. The tur- 
bulence includes the forces from the whitecaps, the shear from the breaking wave. 
internal waves, and Langmur circulation. Thorpe’s equation representing the process 1s 
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Where 1, 1s the rising velocity of the bubbles in msec, N is the number of bubbles per 
unit volume, KX, is the turbulent diffusron coefficient m n7?/ sec , and o 1s the bubble loss 
rate due to diffusion in number of bubbles per second. 

There have been many other techniques developed to measure these background 
bubbles. Johnson and Cooke (1979) and Toba (1961) utihzed a photographic method 
to accomplish this task. Three tmages corresponding to the specular reflections of the 
light provided by three strobes served to idenufy bubbles and as the basis for measure- 
ment. The depth of field was established by the intersection of three dots and a high re- 
solution film allowed a 3X magnification. This method claimed identification of bubbles 
with radii as small as 17 wm. The procedure requires the water to be in a clean, clear 
chute and can identifv bubble radii with only poor resolution. The abundance of debris 
and other matter would interfere with accurate identification of bubbles with a radius 
Sigless than Cyl nnn, solovavev (1976) used a bubble trap in order to concentrate the 
bubbles to allow measurements by photographic means. Breitz and Medwim (1989) de- 
veloped a floating acoustical resonator to determine the number and size of bubbles 
produced from spilling breakers at sea. The Q of the resonator was changed by bubble 
presence. The system allowed measurement of bubble populations of nine radn between 
30 and 270 wm. Though accurate, this method measured bubbles at a depth of 25 cm 
below the surface and getting densities at a shallower depth mav not be possible. 

With these successful results in measuring and understanding the background bub- 
festive ext major area for study is to exanune the origins of the bubbles and the re- 
lationships between these bubble producing actions. 

A large proportion of bubbles in the oceans originate from the whitecaps of break- 
mee vaves. |lolthuijsen and Ilerbers (1986) conducted measurements of wave herght, 
Period, steepness, crest front steepness. horizontal crest asvminetry, and vertical crest 
asymmetry, in an attempt to estimate from these parameters the likelihood that a wave 
will break. Thev concluded that breaking waves cannot be separated from non-breaking 
waves with these factors. Thorpe and Ilall (1983) arrived at several conclusions from 
observations of waves at sea: separate bubble clouds are created by individual waves; the 
breaking crests are seldom more than several meters long: waves generally continue to 
mean as tliey advance so the horizontal size of the bubble clouds increases with wind 
speed. Another mechanism of bubble production involves the small capillary waves that 
€an be observed on that front face of long gravity waves even in the absence of wind. 
These capillary waves can break and pinch off air bubbles in the same manner as the 


breaking gravity waves. 


The role of surface tension in breaking waves is also important. Muller (1972) per- 
formed a series of experiments on lowering the surface tension of the water with chem- 
cals. In these experiments it was found that, as the surface tension decreased, the waves 
got higher before they broke. Gucinski (1986) found that the persistence of visible bub- 
bles at sea ts longer than that in fresh water during similar wind and sea conditions. A 
contributing factor is the presence of organic matter which surrounds the bubble inter- 
faeces. Laboratory eXperiments show an increasing coalescence with increasing 1onic 
SUreig tins 

Since distributions of wind and wind waves can be collected by satellite sensors, the 
development of expressions for related phenomena as functions of these winds and 
waves should be useful. Many studies have been made and great quantities of data have 
been taken of wind stress and whitecap coverage measured from a level Of M"meiens 
Toba and Koga (1968) have shown that the percentage of breaking waves ( @ ). per- 
centage of whitecap coverage ( P ), and the concentration of sea salt particles (eae 


related by the following term: 
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where uw, 1s the surface friction velocity of air in m sec, v 1s the kinematic viscosa 
air in air/ sec, and @ is the spectral peak frequency of the wind waves. |i is mO\geaiee 
possible to measure the tropospheric aerosols over the oceans with satellites. The 
method relates the upwelling visible radiants measured by the satellite to the aerosol 
opucal thickness of the atmosphere. There is a relauionship between the radiance and the 
aerosol optical thickness which varies with parameters such as aerosol size, distribution, 
and refractive index. The correlation between the sources of the aerosols, suchas 
breaking waves and bursting bubbles, and their concentration in the atmosphere ts still 
not understood. 

The use of passive acoustics and modern signal processing allows individual bubbles 
to be exanuned as thev are created bv the breaking wave. Mledwin and Beakyagities 


separated the characteristics of these bubbles into four distinct categories. These 


categories describe the bubble damping, resonance frequency, and oscillation patterns. 
With the method described in this thesis, utuhzing the same equipment as Medwin and 
Beaky, the location of these individual bubbles can be determined from the geometry of 
the hvdrophone array and the difference in the time of arrival of the pressure ficld to the 
individual elements. From this data the exponentially decreasing temporal rate of bubble 
creation from an individual breaking wave can be calculated, as well as the production 
density i bubbles per unit area originating from a single breaking wave. The volume 
of gas encapsulation per unit surface area of a spilling breaker can be computed from 
the bubble density. This could be useful in studies of gas exchange at the ocean surface. 

If all the information and methods described can be integrated and successfully uti- 
lized it will be possible to obtarn a complete perspective on the generation of bubbles 
at sea. The background bubbles are understood and can be measured; the main source 
of bubble creation, parucularly whitecap coverage from breaking waves. can be estI- 
mated; the physical and ehemical makeup of the water can be examined; and now with 
the use of the method described in this thesis it is possible to predict the rate of bubble 
production from the whitecaps. This information, when combined, will further aid in our 


understanding the effect of air’sea teraction on our environment. 


It. FACILITIES AND APPARATUS 


The laboratory facihties include a 57x4x4 foot wave tunnel eniptving =mtom 
10x10x10 foot anechoic tank. A motor-driven reciprocatiig wedge generates the waves 
at the beginning of the tunnel. The waves break approximately three fect after entering 
the tank. The exact position of where the waves will break was not controlled. The tank 
sides and bottoms are lined with Redwood pilings. After spiiling, the waves are absorbed 
by a “beach” consisting of an aluminum shavings Wedge which minimizes the wave re- 
flection at the tank boundary. 

Medwin and Beakyv (1989) conducted tests in the tank which demonstrated at least 
a tS decibel difference between a signal and the strongest reflection from the tank walls 
at frequencies under 5 kI{z. and at least 20 dB at frequencies above 5 kilz. wie 
hydrophone 1s 24 cm from the bubble source. The stgnal to reverberant noise ratio 1s 
better for shallower hvdrophones. 

Two vertical arrays of 2 omnidirectional hydrophones with the shallow hydrophone 
at a depth of 15 cm and the deep hydrophone at a depth of 31 em are situated m the 
approximate area where the majortty of the waves break as depicted m figure {ae 
output srenals are mput mto an Ithaco 1250 voltage amplifier bandpass filter with the 
bandpass frequencies set between 300 [1z and 100 kIIz. The amplifier output splits into 
a Hewlett-Packard 356LA signal analyzer and also into an IBM) AT personal computer 
with an R. C. Electronics “computerscope” analog to digital converter. The signal anaes 
Ivzer and the “eomputerscope’ allow the operator to select triggering level. Ihe Sramma 
analvzer triggers at an input level of 0.837 Volts which. for an amplification of 5000, 
corresponds to an acoustic pressure of I.4 Pascals. 

This “computerscope” analog to digital conversion software package samples the 
data with a sampling period of 1 microsecond with one input, 2 microseconds with two 
miputs, and accordingly, 4 microseconds with four mputs with a }2 bit amplitude resol- 
ution. The sampled signal ts displaved both graplucally and nunterically with curser 
movement. The buffer 1s 64 millrseconds long with the minimum sampling period of 1 
microsecond and can be lengthened by a corresponding increase in the sampling pertod. 


Sections of the buffer can be “expanded” or “compressed” in time and the amplitude of 


the signal displav can also be adjusted. This enables the user to manipulate the graphical 
presentation to obtain a signal that can be accurately analyzed and also to adjust the 


sainpling rate to maxiniuze the signal resolution. 
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Figure 1. Equipment setup. 


IV. EXPERIMENT 


A. BREAKER SPECTRUM 

During World War II] many observations of the ambient noise of the oceans showed 
that “wind noise” is the primary source of the noise in the frequency range between 500 
Ilz and 20 kI1z. The famous Knudsen spectra broke down this large collection of data 
In a manner intended to demonstrate the wind speed dependence of the noise level. Thrs 
spectrum is reproduced from Urick (1983) as figure (2). The possible processes, as given 
by Lrick, which produced this noise were hypothesized to be “crash noises of breaking 
waves . “flow noise” produced by wind moving over the water, “cavitation noise”, and 
“wave generating actions” of the wind. The paper by Medwin and Beakv (1989) proves 
that the catastrophically created, damped bubbles, by themselves, explain the Knudsen 
spectra. The slope of the noise spectrum in this frequency band 1s approximately -5 to 
-6 dB per octave or -17 dB per decade with the level increasing with wind speed. The 
noise in the frequencies outside this band was believed by Urick to origmate from 
sources not directly dependant on wind speed such as shipping below 500 Iz and ther- 
mial noise above 50 KI Iz. 

In order to deternune the norse spectrum of the waves breaking in the tank it was 
first necessary to ensure that the wave generating equipment was not acting as an addi- 
tional noise source. Figure (3) shows the background ambient noise of the tank with no 
equipment in operation plotted with a dashed line. The ambrent noise in the tank with 
the plunger in operation but no waves breaking is plotted with a solid line. Both plots 
were generated on the HP 356!A spectrum analyzer using a ttme average of 30U Input 
stenals, cach of 500 millrseconds duration. The signal level was computed for frequencies 
iperevecn S00 1iz and 50 kFlz. 

The breaker spectrum was also obtained using the TIP 3561A with the input from 
an ommidirectional hydrophone at a depth of 24 cm. The signal was divided into 40 se- 
quential records, each of 20 millisecond duration. Each 20 nullisecond band was Fourier 
transformed to provide the power for each 50 T1z block of the bufler.* The energy in 
these bands was summed only if the energy was 3 dB above the noted ambient noise level 
of the tank with no activity. This was done to ensure that the energy originated onlv 
from breaking waves in the tank. The breaking of a wave can usually be identified bv 


a relatively high amplitude low frequency bubble which served as the triggering source 


for the signal analvzer. The frequencies that were observed were restricted by the spec- 
trum analyzer to between 400 I1z and 20 KHz. A total of six breaking waves was meas- 
ured and the levels averaged. The results are displaved in figure (4). The slope of the 
spectrum 1s similar to the Knudsen spectra at approximately -17 dB per decade. The 
norse level at a frequency of 1 kHz is approximately 82 dB re lpPa?/J/z which 1s some- 
what greater than that of the Knudsen spectra at that frequency at high wind speeds. 
This is due to the laboratory measurements being made only when and where bubbles 
are being produced: a comparable situation at sea would require breaking waves to exist 
everywhere at all times as discussed by Medwin and Beakyv (1989). The low noiseuieuy 
of figure (3) confirmed that the only significant sources that contribute to the sound 
levels above the ambient room noise are breakers produced by wave actions and not 
mechanical equipment. It should be observed that this phenomenon was observed under 
the condition of no wind over the tank. This confirms the suspicion that the actual 
source of the ambient noise are the breaking waves and not the wind. It was the fasci- 
nating fact that the spectrum in the tank is very similar to the expected spectrum in the 
oceans that convinced the author that the data produced by this study may be of future 
use to actual operations performed at sea. The method of using wind speed as a reference 
is useful, however, because the wind speed can be more dependably measured than the 
sea state. An interesting thing to note about the breaker spectrum is that below IR 
the level remains constant to the minimum frequency measured at 400 [1z. 

M.S. Longuet-Ihggins predicted that certam frequencies would be predominant 
from bubbles produced by breaking waves. These frequencies originate from modal dis- 
tortions in the shape of the bubble durmg its creation. The frequenctes were stated in a 
subnutted paper, Longuet-Iliggins (1989), to be at 1.38, 1.69, 2.10. 2.60, 3.32) 4333 ieee 
7.94, T}.4.and 16.8 kflz. These freguencies are designated by tangles on inure (4 je 
laboratory spectrum showed a peak onlv at two of these frequencies, 1.69 and 2.60 KIIz. 
Several of the other frequencies correspond to the position of a minimum) on tie spec 


trum. This does not allow a determination of the significance of this phenomenon. 
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Figure 4. Average noise spectrum from six breaking waves: Triangles mark the 
frequency peaks predicted by Longuet-Higgins. (1989) 


B. BUBBLE IDENTIFICATION 

It is possible to positively identify individual bubbles within the sound originating 
from the breaking of a wave. The equipment utilized is a hydrophone. a preamplifier, 
and an IBNf AT personal computer with an R.C. Electronics Computerscope analog to 
digital converter. The converter samples the signal and allows detailed viewing of the 
desired segments of the captured noise. 

When a wave breaks or “spills” air is entrained, bubbles are formed, and the bubbles 
oscillate to become sources of noise. In every wave observed with the signal processing 
equipment, the moment when the wave breaks can be identified by a large amplitude, 
low frequency bubble, with sufficient energy to be used as the triggering source for the 
signal analyzing equipment. 

This equipment has already been used by Medwin and Beaky (1989) to characterize 
the bubbles from breaking waves by identifving the damping and modulation aspects 
most commonly observed. Four distinct bubble scenarios were classified: spherical bub- 
bles with either a constant decay rate or with two different decay rates, a damped oscil- 
lation bubble appearing to shed a much smaller bubble. a near surface bubble which has 
an increasing amplitude as it moves away from the surface due to an merease igtlie 
dipole axis, and a bubble undergoing amphtude modulation possibly due to another 
bubble in close proximity with a shghtly different frequency. 

The resonance frequency of the bubbles was converted to bubble radius using the 
following relauonship from Clay and Medwin (1977): 
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where fis the resonance frequency of the bubble in Ez, ais the bubble radius in meters, 
y is the ratio of specific heats of the gas in the bubble, P is the hydrostatic pressure in 
Nit, and op 1S ic Ucisit\ Ofte \ ater ines cae. 

The large amplitude initial bubble created by a breaking wave can be seen in figure 
(5). This is the first 120 milliseconds of a S00 millisecond buffer. In this printout the 
upper and lower signals are from a hydrophone at 16 em depth and a hydrophone at 31 
cm depth respectively, The output signals are 1$0 degrees out of phase because ost 
amplifier output phase difference. The confirmation that it 1s indeed a bubble can be 


made from the damping relauonship from Medwin and Beaky (1989) 
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mMliere t 1s the decay time in seconds to | e of its amplitude, fis the resonance frequency 
in Hertz. and 6 is the damping constant. The measured frequency of the bubble in figure 
isso U5 lz and the measured |e decay aime ts 30.9 mulliseconds. From the table on 
page 199 of Clay and Medwin (1977) the damping constant for this frequency should 
be approximately 0.022. If placed in the above equation this damping constant yields a 
Merdecay time of 28.6 milhseconds, a difference from the measured value of 7 percent. 
This is close enough to the theoretical value to ascertain that the bubbles are behaving 
as expected. If the presentation is expanded and viewed with more detail as in figure (6), 
additional later signals from higher frequency, smaller amplitude bubbles can be seen 
superimposed on the low frequency bubble signal. Finally, the presentation can be ex- 
fomced even further, figure (7), to enable the accurate identification of a smaller bubble 
with frequency 16.6 kf fz (196 micron radius) and accurate measurement of the difference 
in time of arrival of the signal to the individual hvdrophones. 

It can be seen in the above figures that the mdividual bubbles can be identified with 
enough resolution to obtam accurate information about bubble frequency, amphtude, 
and relative ume of arrival at the hydrophones. The signal data was recorded onlv ifa 
bubble could be confidently identified above anv noise which may have been received. 
With this equipment accurate identification could be ensured if the signal level was at 
eames decibels above the measured ambient noise level of 37 dB re InPa/i/z at the 
i arepvonesm= sie resulting bubble count ts thereby a conservative measure of the 
nuniber of bubbles and would not include any bubbles with received signal less than the 
hinut of 37 dB. 
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Figure 5. 


04.8 wS/Div 128, 688 uS 


Large triggering bubble from breaking wave: Predominantly large 
amplitude and low frequency bubble (505 Hz). Signal presentation is 120 
ms. Outputs from two hydrophones shown. One output is inverted by the 


amplifier. Sampling period is 15 microseconds. 
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Figure 6. Additional smaller bubble identification: Several bubbles with higher 
frequency can be seen superimposed on the initial bubble signal. This is 
15 ms of the 120 ms plotted in figure 5. Sampling penod is 15 microsec- 


onds. 
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Figure 7. Small bubble identification: Frequency is 16.6 kHz. This is a 1.875 
ms section of the signal shown in figure 5. Sampling period is 15 micro- 


seconds. 


C. BUBBLE LOCATION 

The location of the bubble can be deternuned using the geometry of the hydrophone 
positions and the relationships developed from the difference in time of arrival of the 
sound pressure wave to the hvdrophones. In this procedure it was assumed that the 
bubbles were created and remain approximately at the surface of the water. The argu- 
meit for this assumpuon follows: the bubbles were produced from the action of the 
breaking wave entraining pockets of air. This catastrophic force was the driving mech- 
anism causing the oscillation. It was observed that the bubble oscillations only lasted for 
a short duration, of the order of tens of milliseconds, after the bubble formation. In thus 
short period of action the bubble could not physically travel any signifteant distance. 

A second assumption 1s that the surface of the water remained horizontal, tgnoring 
item ave hereht. | he error due to this assumption will be exanuned shortly. 

When two hydrophones are placed in a vertical array the sound pressure wave from 
Emoupble at tlie suriace will arrive at the hivdropliones at different times. This difference 
in tune of arrival will be a function of the distance the wave must travel through the 
Beater. 

Bubbles produced near the surface of the water act as dipole sources. The dipole 
consists of the bubble and an imaginarv source the same distance above the surface as 
the bubble 1s below. with a distance L separating this image and the bubble. With the 
assumption that kL < 1, the dipole radiation pattern will be defined by the following 


equation from Morse and Ingard (1968): 
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a “s JocD cos(@)(t + an ern 
Where k ts the wave number in radians'meter. R is the range in meters from the bubble 
to the hydrophone, pc ts the acoustic n:mpedance of water in Pa sec.m, D = VL where 
= 4dreU is the volume per time. U ts the radial velocity amplitude. @ is the angle be- 
tween the axts of the dipole and the line adjoining the bubble and the hydrophone, and 
wis the angular frequency of the oscillating bubble in radians second. The magnitude 


of the pressure at hydrophone | can be described as 
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The real part of the pressure at hydrophone | can now be described by the equation: 
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Identically, the real part of the pressure at hydrophone 2 can be written: 
ee ee , . T> ie —! ] 
Pe=i is ecOstil, — kK, tale eee): 
7 z ~ Ak, 


Comparing the phases for the times that the peaks occur, 








Si! 
Wt) vn) fe La at ER, ? = 717 
and 
Wh — KR + tan” (FR) = ne 
and subtracting vields: 
il | —| l 
wii, — 1, ) = ACS — 1) Sea ae oe yeu (ER? =. 


The equation describing the difference in time of arrival can be written as: 
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where R, and R, are the distances from the bubble to the hvdrophones in meters, k 1s the 
wave number in radians meter, c is the speed of sound in meters ‘second, and f 1s the 
resonance frequency of the bubble 1n [lertz. 
The geometry of the hyvdrophones is depicted in figure ($). The following relation- 


ships should be noted: 
eee 
R, =(hy + x°)2 
7 a 
R, = (hy 4+ x") 2 


Where fi, and A, are the depths of the hvdrophones in meters and x 1s the horizontal dis- 
tance on the surface of the water above the hydrophone array in meters. If the position 
on the surface above the hvdrophones is treated as an origin. x 1s in fact the radius of a 
omelet: these equations for K, and K, are placed in the previous equation for time of 
arrival difference, this radius, x, will remain the only unknown. Using an iterative pro- 
gram, the computer will easily deternune the solution. 

The position of the bubble on the surface can be determined by using two vertical 
arrays of hydrophones. The bubble will be located somewhere on the circles surrounding 
Gierarrays. | hese circles will tntersect at two points, in a symmetrical manner about the 
horizontal line separating the arravs as in figure (9). The ambiguity between the two 
Intersections 1s unimportant to the objectives of the experiment. 

The fact that the surface of the water is not constantly horizontal, but in reality 
consists of various wave heights. causes a change m the time of arrival difference be- 
pyeen the hyvdrophones of an arrav. I{ the bubble 1s located o1 a wave crest or trough 


the distance the pressure Wave will have to travel can be described by: 
2, 244 
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Where Af is the vertical difference in bubble position from the horizontal surface of the 
water. This difference in time of arrival caused by the wave amplitude is shown in figure 


(10). It will be noted that the maximum amphtude of the waves in the tank is 


ns 


approximately 3 cm. At this wave height the error in time of arrival is only 5 microsec- 
onds. Because of the fact that the wave herght follows a Gaussian distribution, the wave 
height at the bubble location will likely be on the order of 1 cm and the error mm trme 
would be 2 microseconds with a corresponding error m bubble location of an insignif- 
rcant 2 or 3 millimeters. 

Two vertical arravs, each including a shallow hydrophone at 16 cm depth, and a 
deep hydrophone at 31 cm depth, were placed in the tank in the area where most waves 
break. The perpendicular line between the hydrophone arrays was set at approximately 
a 30 degree angle to the direction of the wave propogation im order to ensure geomeme 
efficiency in the calculation of bubble positions. The buffer length was set to 100 mullt- 
seconds with a corresponding saimpling period of 6 microseconds (spatial resolution, 9 
mm). This small sampling period was considered necessarv in order to measure the time 
of arrival difference as accurately as possible. The output went through the preamplifiers 
set at a gain of SU00 and acted as the trigger to the computerscope analog to digital 
CONVEMer, 

The difference in time of arrival to the two vertical pairs of hydrophones, the reso- 
nance frequency of the bubble, and the time that the bubble began oscillating relative 
to the breaking of the wave, Were recorded in the first 100 milliseconds for each bubble 
created by a breakme wave. This information was used with the iterative compas 
_ program to plot the locations of the individual bubbles. Figure (11) is the result of one 
breaking wave showme the svnunetrical distribution of the bubbles created in the first 
100 nutitseconds alter breaking. The symmetrical locations are caused bv the two inter- 
sections of the bubble radn circles about the axes of the hvdrophones. As lone asutie 
wave does not break close to the hydrophones it 1s not necessary to plot both sets of 
location positions because 101s the relative positions of the bubbles that 1s necessary for 
density calculations. figures (12) and (13) show the bubble locations with only the pos- 
WIMCMPOSILigns TeECcondca, 

The positron plots demonstrate the bubble patches or “hot spots” often noted in 
acoustical measurements at sea. The plots were divided bv a grid consisting of squares 
5 cm ona side. These blocks of 25 square centimeters were used to calculate tieyages 
over Which bubbles were produced. This large area assigned to any bubble counters aia 
compensates for anv error in the calculation of bubble position. Each square with a 
bubble in residence was counted. The average total area of bubble production was 345 
square centimeters per breaking wave with a standard deviation of 42 square centinie- 


ters. It was assumed that ali bubbles created bv a breaker origmate in such a calculated 


area. It was not possible to expand the time of observation in order to confirm this as- 
sumption for later bubbles because the sampling pertod would also have become longer 
and the resolution of the time of arrival difference would not have been sufficient for 
accurate bubble location. Another method was to form a rectangle around the bubbles 
and compute the production area. The average production area calculated with this 
method was 320 square centimeters with a standard deviation of 96 square centimeters. 
The rectangle is marked by a solid line around the bubbles tn figures (12) and (13). The 
two methods compare within 10 percent for the average production areas. 

Another result of the bubble position plots was the abihtv to label the order in which 
the bubbles were created by the breaking wave. Figures (12) and (13) also show these 
relative positions of the bubbles from two waves each marked with the order in which 
they were produced. The order of the positions does not demonstrate that the location 
of the bubble depends on the relative tine of its creation with respect to the breaking 
@ieuie wave. Ihe jirst bubble created was not necessarily the closest to the position 
where the wave breaks. Also the last bubble created was not necessarily located the far- 
thest along the direction of wave propagation from the position where the wave began 
meaning. Iins fact suggests that the majority of the bubbles created were located within 
the area calculated even though this area was calculated from bubbles originating in the 
first 1U0 mulliseconds of breaking. 

Bubbles observed within the bubble location plots were ocasstonally positioned in 
roughly linear patterns not necessarily perpendicular to the direction of wave 
propogation. The bubbles within these patterns all originate consecutively in short 
Pidows of tine. It is believed that these bubbles are created by and along the small 
capillary waves which sometimes are seen superimposed on the face of a main wave. This 


evidence of capillary waves is examined in Appendix A. 
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Figure 8. Bubble location geometry. 
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Figure 9. Bubble location technique: Intersection of the two circles of radi x _ 
and x, around the hydrophone axes marks the possible location of the 
bubbles. 
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Figure 10. Timing error due to wave height: Error in time of armval caused by 
the additional distance traveled by the noise of a bubble on the surface 


of the water due to local wave height above mean surface. 
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Figure 11. Distribution of bubble locations: Distribution is symmetréc about the 
line joining the two vertical arrays. Bubbles are from the first 100 ms 
after the breaking of the wave. 
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Figure 12. Bubble location plot: Bubble locations marked with the order of bubble 
creation. Same breaker as figure 11. Bubbles are from the first 100 ms 
after the breaking of the wave. 
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Figure 13. Bubble location plot: Bubble locations marked with the order of bubble 
creation. Bubbles are from the first 100 ms after the breaking of the wave. 


D. RATE OF BUBBLE PRODUCTION 

It was observed from the use of the Hewlett-Packard 3561A Spectrum Analyzer that 
noise is produced from a breaking wave for a period of 2 to 3 seconds aiter thiemijms 
sound of a breaking wave is heard. The digital processing equipment used to identify, 
locate, and count bubbles has a relatively small memorv (64k bites). This does not per- 
mut the sound to be captured for a long enough period of time to record all bubbles 
produced and simultaneously have a sampling period small enogh to be of use in the 
identification of bubble resonance frequencies. Therefore it was thus necessary to deter- 
mune the rate at which bubbles are produced over time in order to determinewamees 
pression defining the bubble production rate relative to the time at which the wave 
breaks. The measurements could then be made for shorter periods of time and then ex- 
trapolated to determine the total number of bubbles produced. 

Two omnidirectional hyvdrophones Were placed in a vertical array in the tank at 
depths of 12 and 24 centimeters. The buffer length for two hydrophones was expanded 
to SOO milliseconds by using a sampling period of 15 mucroseconds. Tlirough exper- 
Imentation and consideration of the Nyquist relationship 1t was determined that this was 
the maximum sampling period that would allow dependable identification of bubbles 
with resonance frequencies less than 33 kHz. The trigger level was adjusted 1Gyamigaas! 
sufficiently high to ensure the sampling equipment would begin capture only when a 
wave was breaking almost directly above the hydrophones. This was done to guarantee 
the observation of bubbles with a sinall amphtude which may not have been seen from 
a longer distance and at a greater angle from the dipole. The short range from the 
hvdrophones can be confirmed using the method previously described. The input signal 
passed through the preamplifier with the gain set on 5000. A total of 10 breaking waves 
were Observed. 

The average number of bubbles of all radit produced per breaker, within time incre- 
ments of 10 nulliseconds, 1s plotted against the time after the breaking of the waveum 
figure (14). Additionally. the production rates for bubbles with bubble radu broken down 
into bands of greater than 2.2 mm (resonance frequency less than 1500 HZ), between 22 
mum (1500 Hz) and 0.16 mm (20 kJiz). and less than 0.16 mm (greater than 20 KHZ), are 
depicted in figures (15). (16). and (17) respectively. These bands are chosen to accentuate 
the bubble production behavior in the Knudsen sea noise spectrum frequency regime 
senerilly taken to be betwee S00 112 ies. 

The production spectrum of the two radi bands, 0.16 to 2.2 mm and smaller than 


0.16 mm as seen in figures (16) and (17). suggests that bubbles in these bands are 


produced with an exponential rate of decrease. The equations which approximate the 


decay rates for these two bands are: 


aN = 0.3] e707 
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and 
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respectively. This 1s probably also true for the largest bubbles but not enough data ts 
available to show it. 
The curve which best approximates the exponentially declining total production rate 


of figure (14) ts: 





where dN ‘dt 1s the number of bubbles produced per mulltsecond and t ts trme after the 
breaking of the wave in milliseconds. The expression stated above can be integrated to 
describe the bubble production rate out to specified trmes. When the integration was 
performed to infinity it was found that less than 0.1 percent of the bubbles would origi- 
nate after the first second. About 97 percent of the bubbles were formed in the first S00 


maliseconds alter the wave breaks. 
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Figure 14. Average number of bubbles produced per laboratory breaker: (all 


sizes). The bubbles are counted within 10 ms increments. 
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Figure 15. Average number of large bubbles per laboratory breaker: (radius > 
2.2 mm). The bubbles are counted within 10 ms increments. 
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Figure 16. Average number of mid-size bubbles per laboratory breaker: (0.16 


mm radius<¢ 2.2 mm). The bubbles are counted within 10 ms increments. 
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Figure 17. Average number of small bubbles per laborat. ~ breaker: (radius < 
0.16 mm). The bubbles are counted within 10 n_ .ncrements. 
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E. BUBBLE PRODUCTION DENSITY 

The rate at which bubbles will be produced is related to the foam coverage. With 
this knowledge researchers will be better able to understand the phenomena which affect 
acoustics, Weather. aerosol production, and any other sea'surface interactions. The two 
major requirements for density calculations have already been covered: the bubble cre- 
ation rate relative to the start of the breaking of the wave, and the area over which these 
bubbles are being produced. 

All bubbles in figures (15) - (17) with equal radi) were summed and the average 
number of bubbles per breaker with that radius was computed by dividing by the num- 
ber of waves observed. In an earlier experiment itt was found that only 3 percentage 
bubbles originated in the second 500 milliseconds after the breaking of the wave. Because 
the measurements were again taken over a period of 500 nulhseconds the numberor 
bubbles observed was thus increased by 3 pereent to give the number of bubbles in the 
first second. This number 1s now defined as the “total” number of bubbles. In an addi- 
tional earlier experiment it was found that the average area of bubble production 
spanned 345 square centimeters or equivalently 0.0345 square meters. The “total” 
number of bubbles was divided by 0.0345 to give the resulting average bubble production 
Jensity in bubbles per square meter for each bubble radius observed. Thus the bubble 
production will be expressed in bubbles per square meter and assumes negligible num- 
bers of bubbles are formed bevond the first second after the wave breaks. 

Bubbles were observed with finite numbers of diflerent size 1adu, If an mnie 
number of waves were measured, eventually bubbles would be observed having every 
possible radius up to the maximum gas volunie that can physically be encapsulated. To 
account for this mutation in the data recorded, the span m bubble radius in microns 
was calculated to include half the bubble radius size between each of the subsequent 
smaller and larger bubble observations. The number of bubbles per square meter was 
divided by this radius span to give the bubble density in bubbles per square meter per 
mucron of bubble radius increment. 

The bubble production density 1s plotted on senulog paper with the radius on a 
logarithmic scale in figure (18). and both bubble radius and density on a logarithmic 
scale in figure (19). It can be seen that the most conunonly produced bubble was 0.15 
mum in radius produced at the approximate density of 5.9 bubbles per square meter per 
nucron radius increment. The smallest bubble radius measured was 0.048 mm which has 


a corresponding resonance frequency of 66 klIz. The largest bubble was of radius 7.4 


mm with a corresponding resonance frequency of 440 iz. It will be assumed that still 
larger bubbles would not have been found even if more breakers had been observed. 
Once the production density of bubbles for each radius has been calculated the vol- 
ume of gas encapsulated from each known bubble radius can be computed. Figure (20) 
is a summary of cubic centimeters of gas volume encapsulated per square meter of 
breaking wave for bubbles of various radi. Although it was seen in figures (18) and (19) 
that the majority of bubbles produced were less than 0.5 mm in radius, the contribution 
to total gas entrainment from these small bubbles 1s minimal. Assunung that there are 
no bubbles larger than the largest measured in 10 breakers which had a radius of 7.4 
mm, the individual contributions from each bubble radius can be summed to give a 
mele Tate Of eas “encapsulation which averages 25 cubic centimeters of gas 


encapsulated per square meter of breaking wave in the laboratory. More than half of 


tus total volume came from bubbles with radius greater than 3 nullimeters. 
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Figure 18. Average surface density of bubbles produced. 
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Figure 20. Volume of air encapsulation per surface area. 
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V. COMPARISON WITH OTHER DETERMINATIONS 


We know of only one other attempt to measure the bubble production in a spilling 
breaker, that was by Toba (1961). Contrary to our breakers being produced by an os- 
eilanne wecdse wath no wind, }oba used wind of different speeds and fetch in a flume. 
Toba’s method considered that an equilibrium existed between the bubbles entering from 
periodic breaking at the water surface and those rising to the surface. Fis photographs 
were made of the volume density of bubbles, not the surface production density of the 
femmes Which 1s OUr goul, Ile used the time between breakers to infer the bubble 


encapsulation rate using the relationship 





x 
| 100 ae 


~~ 


100 
centage of crests actually entraming bubbles. © 1s the total number of bubbles of a par- 





where fis the bubble encapsulation rate in bubbles ‘square cm second, 1s the per- 
ticular size m a Water columin of unit base area, and T 1s the period of the maim surface 
eves lhe shape of Toba s encapsulation spectrum, heures (2)Pyand (22). was very sim- 
Herero tliat in this experiment with a sharp peak at a single bubble radius. However the 
peak of Toba’s bubble production was at approxmmately 0.38 mm. in radius at a rate of 
meproatnatcly 0.2 bubbles per square centineter per second for a wind speed 12.1 m’‘sec. 
[lis photographic measurements mcluded 6 different sizes of bubbles with radu from 0.15 
mim to 6.0 mm. Our method observed 53 bubble radi: from 0.048 mm to 7.4 mm. Flis 
smallest radius bubble was our most common bubble. The production rate inferred by 
Momator 12.1 msec wind was of the same order as the rate measured directly acous- 
tically. 

Micre diac been extensive nicasuremients of the foam coverage at sca usmg aerial 
photography and satellite maging. The foain coverage is 2 measurement of residual 
bubbles including those convected upward from below and lacking those that have been 
Semvected from the surface to lower depths or that have burst at the surface. Thus the 
areca of coverage can not be used directly with the densities measured m this experiment 
to calculate the bubble populations beneath the ocean surface. Further experrmentation 
is required in order to understand exactly what ts being measured from the foam of 


breaking waves. 
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Figure 21. Rate of formation of air bubbles in wind waves (Toba). 
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Figure 22. Rate of formation of air bubbles in wind waves (Toba): plotted per 
equal log (bubble size) interval. 
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VI. CONCLUSIONS 


The method described throughout this thesis achieved its goal of obtaining the 
bubble production density from plunger generated spilling breakers. The technique used 
was to passively observe the amphtude and frequency of the randomly pulsed sound 
emitted by the breakers. The results give more details of the characteristics of bubble 
production than previous experimental techniques used for obtaining surface bubble 
densities. 

Bubbles were positively identified throughout a radius range of 0.048 to 7.4 mulli- 
meters. Most of the bubbles originate in the first 500 nulliseconds after the breaking of 
the wave though sound from bubbles has been recorded for a period of 2 to 3 seconds 
after breaking. The bubbles were produced with an exponentially decreasing rate fol- 
lowing a decay constant of -0.007t with time measured in milliseconds after the breaking 
of the wave. The average area of bubble production from breaking waves in the tank 
was 345 square centineters. The most commonly produced bubble had a radius of 0.15 
nullimeters and was entrained at a density of 5.9 bubbles per square meter per micron 
of bubble radius increment. This bubble radius is smaller than the minimum bubble ra- 
dius measured in other expernnents using photograpliuc techniques to obtain bubble 
volume density. The greater resolution using acoustics instead of photography has al- 
towed 10 times the resolution in bubble radu idenufication for bubbles with radu be- 
tween 0.25 and 2.0 mullineters. The average encapsulauon of gas was 23 seme 
cenuumeters per square meter of breaking wave. 

Lsing the Knudsen sound radiated from the breaking waves as our sole criterion, 
Whitecaps from spilling ocean waves have htthe physical difference from the breaking 
waves produced in the laboratory. If this is the case, the measurements reported wire 
are useable for the ocean. The use of this passive acoustical technique i aciiaieae 
measurements may also be possible. If these laboratory data are confirmed by actual 
measurements of surface bubble production density at sea, and after further studies of 
bubble dissipation in a turbulent ocean, satellite measurements of the foam coverage 


from whitecaps could vield an estimate of the rate of bubble production in the ocean. 


44 


APPENDIX AxwweEVIDENCE OF CAPILLARY WAVES 


Upon examination of the bubble location plots it was discovered that in some of the 
plots the bubbles appeared to be located along lines perpedicular to the Wave propa- 
gation direction (frg 12) and in other cases the lines are not perpendicular to the direction 
of wave propagation. These lines appeared with varying angles ranging from 30 to 90 
degrees awav from the propagation direction. It was hypothesized that these mav have 
been capillary waves forming on the face of the main wave. A purely arbitrary reference 
trme of 10 mulltseconds was chosen as the time within which all of the bubbles on the line 
must have origmated. Figure (23) is an example of a capillary line identification. The 
bubbles were captured with a sampling period of 6 microseconds during the first 100 
mulhseconds after the wave broke. The order of bubble formation ts marked on the fig- 
ure and the following are the times of the beginning of the individual bubble oscillations 


after the breaking of the wave: 


Bubble nuniber Origin tume (ms) 
| AS 
2 5.88 
3 10,56 
4 lee 
5 1244 
6 18.74 
] 20.54 
S Ne OM, 
9 Otello 
10 BSE00 
Bl 42.37 
IZ 42.49 
13 49.03 
14 Oo. 
IES 69.3 
16 les 
1/7 TSO) 


{. 
‘sr 


1$ 80.45 
19 81.83 
20 93.94 


The capillary lines on figure (23) are idenufied by bubbles 6 through 8 and bv bubbles 


17 through 19. Figures (24) and (25) are photographs of breaking waves in which the 


capillary waves mav be clearly observed. 


+6 


60.0 


hh SS 

eae SaaS Sasa 7 

Re ae ae | fe 
cae eee | 


[S LL 
Ze ae a = 
ee ae a a 








40.0 





y position (cm) 


3 | Sea es hse eee eee 
Cc 
ON) 
| 





BEER EESEE REESE REECE 


' 60. 0 -40.0 =20-9 0.0 20.0 40.0 60.0 


X position (cm) 





Figure 23. Possible capillary wave positions by bubble location: Lines are from 
bubbles 6-8 and 17-19. 
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Figure 24. Photograph of capillary waves. 
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Figure 25. Photograph of capillary waves. 
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APPENDIX B. VOLTAGE ANIPUIFIERS 


When recerving srgnals from many hydrophones simultaneously, it was necessary to 
develop a substitute for the lmmited number of Ithace 1201 Preamplifiers. A series of 
preamplifiers was built using two 356 operational amplifiers mn series with the necessarv 
feedback elements. Figure (26) 1s the diagram of the amplifier circuit. Figure (27) shows 


the operational characteristics of the amplifier. 
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Figure 26. Amplifier circuit diagram. 
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Figure 27. Operational characteristics of the amplifier. 
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